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Abstract: The time-integrated current of the TASEP has non-Gaussian fluctuations of order 
t 1 / 3 . The recently discovered connection to random matrices and the Painleve II Riemann- 
Hilbert problem provides a technique through which we obtain the probability distribution of 
the current fluctuations, in particular their dependence on initial conditions, and the stationary 
two-point function. Some open problems are explained. 



1 Introduction 



The asymmetric simple exclusion process (ASEP) has become a paradigm for nonreversible 
stochastic particle systems. In this note we will consider the particular case of one dimension 
and right jumps only (TASEP). The TASEP has state space {0, 1} Z whose elements, rj, are 
particle configurations with rjj = 0, 1 for j G Z. rjj = 1 means a particle at site j and rjj = 
means no particle, resp. a hole, at site j. The stochastic dynamics is governed by the following 
jump rule: independently each particle jumps with rate 1, i.e. after an exponential waiting time 
with mean 1, to the right neighboring site, provided it is empty. Our rule defines the Markov 
jump process r] t , t > 0. r] is the starting configuration which may be random. 

For the purpose of the introduction let us assume that rj is distributed according to the 
Bernoulli measure /X1/2, i-e. rjoj, j G Z, are independent with P(r/oj = 1) = \- Since fii/ 2 is 
a stationary measure for the TASEP, the process rjtj,t 6l,j G Z, is stationary in space and 
time. We are interested in 

N t = number of particles which have crossed the bond (0, 1) up to time t. 

Clearly, 4fN t is the current across the origin, hence our title. By stationarity, 

E(N t ) = \t (1.1) 

and the real issue are the fluctuations of 

N t -\t. (1.2) 

To convince the reader that the current fluctuations have something interesting to offer, we 
first study the large deviations of N t . For the lower tail we consider ¥(N t < at) with a < |. To 
reduce the current it suffices to let a single particle move more slowly. The other particles will 
then pile up behind. Therefore we expect 

P(A"t < at) ~ e - g - (a)t , a < \ , (1.3) 

for large t with g_ > and <7_(|) = 0. On the other hand for the upper tail, F(N t > at), a > |, 
the current must be increased, which requires order t particles to jump faster. Therefore 

P(At > at) ~ e - 9+{a)t2 , a > | , (1.4) 

for large t with g + > and g+{\) = 0. Since the large deviations have different order of 
magnitude above and below the mean, we must be outside the domain of the central limit 
theorem. 

A more detailed analysis shows that g~(a) ~ c_|a — ^| 3 ^ 2 and g + (a) ~ c + \a — || 3 for a 
close to \ H], 0, also || £| for the PNG model. Extrapolating beyond the validity of the large 
deviation result, we conclude that for the lower tail 

W{N t -\t~x) ~ e - c - |x/ * 1/3|3/2 forx<-l (1.5) 
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and for the upper tail 

F(N t - \t ~ x) ~ e - c+|x/il/313 forx»l. (1.6) 
On this basis the fluctuations should be on the scale t 1 ^ 3 and 

iVt-|t = t 1/3 e (1-7) 

for large t. £ is a nondegenerate random variable with a distribution whose lower tail is 
exp(— c_|x| 3//2 ) and whose upper tail is exp(— c + |x| 3 ). 

The goal of our note is to explain how (|1.7Q and related quantities like the stationary 
two-point function, E(r/o,o^f,j) — f, can be mapped to a last passage percolation problem with 
boundary conditions and possibly defect lines. The asymptotic analysis of such last passage 
percolations has been carried out by Baik and Rains [f| and we will make contact with their 
work towards the end. 



2 Last passage percolation with boundary conditions 

We generalize the set up of the introduction by taking instead of pi/2 the Bernoulli Pp_, p+ 
as starting measure, according to which the r/oj's are independent with P(^o,j — 1) = P- for 
j < and P(?7o,j = 1) = p+ for j > 1. p:y 2 is the special case with p_ = | = p + . Liggett || 
determines the law of T]t as t — > oo. We will come back to his result. Johansson maps the 
special case p_ = 1, p + = to a last passage percolation. We extend his result to arbitrary 
0<p_,p+<l. 

Let us first assign to an ASEP configuration rjtj the height function 

( m + Y,i=i(i-2vt,i), 3>1, 

ht(j)=< 2N t , j = 0, (2.1) 

[ m - El i+1 (i - 2^,0 , j<-i. 

By fiat /io(0) = 0. h t (j) is even at even and odd at odd sites. We will establish that the 
distribution of the heights can be obtained from a last passage percolation problem on the 
positive quadrant NxN with suitable boundary conditions. 

To each site G Nq we associate a random variable w(i,j). Let ( + be geometric with 
parameter 1 — p + , P(C+ = n) = p + (l — p+) n , and independently let £_ be geometric with 
parameter p_, P(C- = n) = (1 — p_)p", n = 0, 1, The {u>(i, j)| (i, j) G Nq} are independent 
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with distributions: 



w(i,j) is exponential with mean 1, i,j>l, 

w{0,0) = 0, 

w(j,0) = 0, 0<j<C+, 

w(j, 0) is exponential with mean (1 — j > C+ , 

w(0,j) = 0, < i < C- , 

w(0,j) is exponential with mean pZ 1 , j > C- ■ (2-2) 

Let Q m , n be the set of all up/right paths on Nq starting at (0, 0) with end point (to, n). To 
each path oj e fi m , n we associate the passage time 

TH= £ (2.3) 

(iJ)6o; 

Then the last passage time, from (0,0) to (to, n), is given by 

G{m,n) = max T(u) . (2.4) 

In the more physical parlance, oo is a directed polymer with end points (0,0) and (m,n). 
Each site carries the energy —w(i,j). The energy of the polymer ou, —T(u), is the sum of the 
site energies along the polymer. Thus —G(m, n) is the minimal energy, or ground state energy, 
of the directed polymer. 

By construction G(m,n) is nondecreasing in both arguments. Therefore the level sets of G 
define a height function which we denote by h t (j). More precisely let A t = {(to, n)\ m,n> 1, 
G(m,n) < t} as a random set. A t is bordered by h t according to A t = {{m,n) \ 2 < m + n < 
h t (m - n)}. 

Theorem 2.1 In the sense of joint distributions we have 

ht(j) = ht(j) for \j\<h t (j). (2.5) 

Proof. The idea is to rotate the TASEP height configuration by — 7r/4 relative to the origin and 
to identify the resulting height differences to the right as a zero range process. Similarly we 
rotate the TASEP height configuration by 7r/4 relative to the origin and identify the resulting 
height differences to the left as a zero range process. We start with the former. 

We define the right dynamics by setting rj ti o = 1 for all t > 0. The rj t j, j > 1, follow the 
TASEP rule. For given rj t ,j let ( t ,j, j > 1, be the interparticle distances to the right. Thus 
if r) tj i = 1, then ( tjl = 0, since the origin is always occupied. Since the starting measure is 
Bernoulli p + for j > 1, at time t = 0, the {Co,j}i>i are independently geometrically distributed 
with parameter 1 — p + . ( t j, j > 1, is a zero range process where the only allowed transitions 
are ((t,j,(t,j+i) ~^ (Ct,j + 1, C*,j'+i — 1) f° r 3 = 1)2,- •• and they occur with rate 1 provided 
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Ctj+i > 0. Ct,i increases in units of 1 and w(j, 0), j = 1, 2, . . . , are the successive waiting times 
of £t,i- At t — 0, Co,i is geometrically distributed as £+■ Thus u>(j, 0) = for 1 < j < ( + . The 
zero range process £tjj j > 2, is stationary. By Burke's theorem the current from 2 to 1 is 
Poisson with intensity 1 — p + . Thus u>(j, 0) are independent exponentials for j > ( + . 

For the left boundary we only have to interchange particles and holes. Holes jump to the 
left only. The left hole dynamics is defined by setting ?7 tjl = for all t > 0. The holes 1 — rjtj, 
j < 0, follow the TASEP rule. ( t ,j, j < 0, are the interhole distances. ( t ,o increases in units of 
1 and w(0, j), j = 1, 2, are the successive waiting times of Ct,o- At t — 0, Co,o is geometrically 
distributed as ( and t \— > ( ti0 is a Poisson process with jump rate p_. 

To build up the random set A t we start with Nq such that (0, j), < j < £+, and (j, 0), < 
j < C- are occupied. For completeness it is useful to regard also the sites (j, — 1), (— l,j), 
j = 0, 1, ... , as occupied. All remaining sites are empty. The set of occupied sites defines 
Aq. The site (m,n) is filled after the random waiting time w(m,n) starting from the instant 
of time when both the left neighbor at (m — l,n) and the lower neighbor at (m,n — 1) are 
filled. A t is the set of occupied sites at time t. Its boundary, h t , is defined in the coordinate 
system rotated by tt/4, to say the coordinates (m,n) are transformed to the new coordinates 
(j, h) through j = m — n, h = m + n. The distribution of h t agrees with h t inside the cone 
C = {(j,h)\ h > \j\}: by construction the events at the boundary of C have the probability 
induced by the TASEP dynamics. In the interior of C the random variables w(i,j) are the 
waiting times for the particle jumps of the TASEP. □ 

Returning to the problem of the introduction, we have to set p_ = | = p + . We define 
t i— > N t as the inverse function to n h G(n,n). By Theorem |2.1| N t = N t in distribution. 
Therefore, if one can control the statistics of the last passage time G(n, n) for large n, one 
knows the statistics of the current across the origin for large t. More specifically, since Nt and 
G(n, n) are linear in average, their fluctuations have the same asymptotic distribution up to a 
linear change in scale. 

3 Defect lines 

There are other cases of interest which can be mapped to a last passage problem. As one 
example we choose the (deterministic) initial configuration as r] j — |(l + (— 1) J '). As before, 
we want to determine the statistics of N t — \t for large t. The last passage refers now to the 
upper triangle A n = {(i, j)\ 1 < i, j < n,i + j > n + 1}. Let Q n be the set of all up/right paths 
starting at the anti-diagonal {(i, j) \ i + j = n + 1} and ending at (n,n). Then 

G pl {n) = maxT(cj). (3.1) 

G pl (n) is the point to line last passage time, in contrast to G(m, n) which is the point to point 
last passage time. As before n \— > G pl (2n) is the inverse function to t i — ^ N t . Odd arguments 
would correspond to the initial condition r] j — |(l — (— 
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One can force the alternating initial condition into the scheme of Section 2. We set p_ = 1, 
p + = 0, and reflect the w(i,j) at the anti-diagonal, i.e. w(i,j) = w{n + 1 — j, n + 1 — i) for 
i + j ^ n + 1. On the anti-diagonal we set w(i,j) = 0. Then in distribution 2G pl (n) = G(n, n), 



where G(n,n) is defined as in ( |2.1|) . In this scheme we have point to point last passage with a 
passage time distribution symmetric relative to the ant i- diagonal. 

A second example is the semi-infinite TASEP with a source at the origin. At site 1 we insert 
a particle with rate a, a > 0, provided site 1 is empty. All other jumps are governed by the 
TASEP rule. The initial measure is Bernoulli p. Then the last passage is restricted to the lower 
triangle {(m,n)\n < m}. The right boundary of the last passage percolation is unchanged. 
The w(i,j), j < i, are exponential with mean 1 and the w(j,j) are exponential with mean a. 
The diagonal is a defect line. To make the formal analogy even closer to the previous cases we 
could copy the w(i,j) in the lower triangle to the one in the upper triangle, which does not 
change the statistics of the last passage time from (0,0) to (n, n). 

A widely studied case is the slow bond problem || As initial measure we take p_ = 1, 
p + = 0. We assume that the jumps through the bond (0, 1) occur with rate r. This corresponds 
again to a defect line along the diagonal, where w(j,j) is exponential with rate r. In contrast 
to the semi-infinite system the w(i, j)'s for i < j and for i > j are now independent. One would 
like to know the average current for t — > oo, joo(^) = lim^oo N t /t. Clearly, joo(^) = r for small 
r and joo( r ) — \ for r > 1. The critical rate r c is the smallest with joo( r c) — \- One knows that 
\ < r c < 1 and conjectures r c = 1. In terms of the directed polymer, at r c there is a depinning 
transition. For small r the directed polymer with end points (1,1) and (n, n) will stay order 1 
close to the diagonal. As r is increased, the size of excursions away from the diagonal increases. 
For r > r c the transverse fluctuations diverge as n 2 / 3 | 1Q| . The directed polymer depins from 
the diagonal. 

4 The two-point function 

The TASEP with Bernoulli /i p (p + = p = pJ) as starting measure is stationary in space-time. 
From a statistical mechanics point of view the central quantity is the two-point function 

S(j,t) =E P (77^770,0) -P 2 . (4.1) 

We list a few properties: 

Proposition 4.1 We have S(j,t) > 0, £\ S(j, t) = p(l - p), £. jS(j, t) = p(l - p)(l - 2p)t, 
and 

8S(j,t) = A(E p (h t (j) 2 ) - E p (h t (j)) 2 ) (4.2) 
with the discrete Laplacian (A/)j = fj + \ — 2fj + and ht(j) defined in ($.\ ). 



6 



Proof: S(J,t)/p(l — p) is the transition probability for a second class particle starting at the 
origin JTTJ, which implies the first two assertions. For the third one we note that 

4 J2jS(j,t) = J2jE p ((r]^ tfi - p)(rj ,j-i(l - Vo,j) ~ Vo,j(l ~ Vo,j+i))) 



dt 

3 3 

= ^2 E p{(V-t,j - P)?7o,o(l - r) 0t i)) = ^E p ((r? j - p)r/ ,o(l ~Vo,i)) 

3 3 

= p(l-p)(l-2p), (4.3) 

where we used conservation of the number of particles. 
To see the fourth identity we compute 

AK p (h t (j) 2 ) = E p (h t (j + l) 2 -2h t (j) 2 + h t (j-l) 2 ) 

= -8M p (N t {r] t!j+1 - THj)) + 2 - 4E p ((l - 2r) 0>j ){r) 0>j+1 - r) 0>j )) , (4.4) 

using (|2.1| ) and stationarity. By translation invariance 

M N t(Vw - Vt,j)) = M( N t~ - N t)Vt, 3 ) , (4.5) 

where iV t ~ is the number of particles jumping through the bond (—1, 0) up to time t. The 
conservation law ensures that 



N t - + ri 0fi = N t + r] t ,o . (4.6) 

Inserting into ( |4.5| ) yields 

E p (N t (rj tJ+1 - rj t j)) = E p ((rj ti0 - Vo,o)Vt,j) (4.7) 

and taking into account that E p (/i t (j)) = 2p(l — p)t + (1 — 2p)j we obtain ( |4.2| ). □ 

(p(l — p)) -1 ^^, t) is a probability distribution with mean (1 — 2p)t One expects [0 that 
its variance grows as t 4//3 which suggests the scaling form 

S(j,t) = p(l - p)(4(p(l - p)) ll H^y l \ g "{(j - (1 - 2p)t)(4(p(l - p)) 1 "^)- 1 ) . (4.8) 

The link to the last passage percolation is provided by ( ^4.2| ) and Theorem which tell us 
that the two-point function is given through the second moment of the last passage time at 
boundary conditions p_ = p, p + = p. Because of the subtracted mean in (|4.8| ) we consider the 
passage from (0,0) to (m,n) for large n with fixed slope 6(p) = n/m = (1 — p(l — p) — (1 — 
2p))/ (1 — p(l — p) + (1 — 2p)). The second moment of this last passage time is proportional to 
g(0) with g the scaling function in (|4.8|) . To have the full function g we have to determine the 
last passage time with an end point of order n 2//3 away from (n/9(p), n). 
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5 The Baik and Rains analysis 



Recently Baik and Rains [0 studied last passage percolation with the above described boundary 
conditions. Unfortunately the case of an exponential distribution is not yet accessible to their 
techniques. Only the result of Johansson is available which in our notation corresponds to 
p_ = 1, p + = 0. For this case he proves that 

hmr 1 / 3 2 4 / 3 (iV t -it) = -^ GUE (5.1) 

in distribution with P(£gue < %) = -Pgue(^)- Fgue is expressed in terms of a solution to the 
Painleve II equation. We will consistently use the notation of || for the various distribution 
functions and will not repeat their definitions here. 

Baik and Rains study the geometric and Poisson last passage percolation. 

(i) geometric. The exponentially distributed random variables w(i,j) are replaced by w q (i,j), 
< q < 1. w q (i,j), i,j > 1, has geometric distribution with parameter q, w q (i,Q), i > 1, has 
geometric distribution with parameter ct+^fq, and w q (0,j), j > 1, has geometric distribution 
with parameter ct-^fq. The random variables {w q (i,j)\i,j > 0} are independent. Baik and 
Rains study the distribution of the passage time G q (n,n) for large n in dependence on q, a + , 
a_, where they allow being close to 1. Specifically they set 

a± = 1 - n~ 1/3 w± (5.2) 

and establish how the asymptotic distribution of G q (n,n) depends on w±. The particular case 



Q!_ = = a + is proved by Johansson ||, who can allow for an arbitrary end point. In |13 
last passage percolation is investigated for a_ = = a + and the w q (i, j)'s satisfying various 
reflection symmetries including the cases of interest here, reflection symmetric relative to the 
diagonal resp. relative to the anti-diagonal. The large n distribution of G q (n,n) is established. 

Following the scheme already explained we see that the geometric case can be interpreted as 
a discrete time TASEP with waiting times w q (i,j) + l. In a single time step every particle which 
has an empty site ahead can jump to the right. They do so independently with probability 
1 — q. Thus the initial step, left density p_ and right density p + , is covered by || and the 
periodic initial configuration, ... 0101 ... , by jnj. To handle the two-point function, a small 
change in the boundary density as in (|5.2j ) should be translated to a small change in the angle 
of the end point for the directed polymer away from the diagonal. This seems to require a 
nontrivial generalization of the known results. 

In the limit q — > 1 the random variables (1 — q)w q (i,j) converge to w(i,j) in distribution. 
Therefore the continuous time TASEP, time t, is recovered from the discrete time TASEP at 
integer times [£/ (1 — q)] as q — > 1. 

(ii) Poisson. We consider the geometric last passage percolation and interpret w q (i,j) as the 
number of points in the square with corners ^Jq{i ± \, j ± |). Then in the limit q — > such 



S 



that (u, v) = yfq (m, n) is fixed we obtain a Poisson point process which has unit density in 
the rectangle [0, u) x [0, v], line density a + on [0, u) x {0}, and line density a_ on {0} x [0,v]. 
An allowed path, u, connects continuously (0,0) to (u,v), and is piecewise linear. Each linear 
segment has a slope 9 such that < 9 < oo and connects two Poisson points, except for the 
first and last piece. T(u) is the number of Poisson point which are transversed by u and, as 
before, 

GWssonK v ) = maxT(u;) . (5.3) 



For a_ = = a + the asymptotic distribution of Gp i SSOn (v, v) is established in |fL4| . Baik and 
Rains H extend to arbitrary and allow for small deviations of the boundary densities 

as in ( |5.2| ) with n -1 / 3 replaced by v _1 / 3 . Very recently they investigate the case of reflection 



symmetry relative to the diagonal, in particular a_ = a> + , and an extra line density along the 



diagonal fllS] . 



The particle model behind the Poisson last passage percolation is the PNG model [fL|], [L7|, [Lq| . 
It consists of point particles with velocity ±1. They annihilate each other at a collision and 
are created in ± pairs with rate 1. If p+{x, t) is the density of right movers, and p-(x, t) of left 
movers, then the field <fi{x,t) = p + (x,t) — p^(x,t) is locally conserved and corresponds to r\ t j 
of the TASEP. More precisely the corresponding height function 

h(x,t)=N t - [ dy<j>(y,t), (5.4) 
Jo 

where as before N t is the time integrated current at the origin, N t = J * <is(p + (0, s) + p-(0, s)). 
The height is related to the last passage time through G(u,v) = h(x,t) in distribution with 
u = (t + x)/y/2, v = (t- x)/V2. 

In contrast to the geometric and exponential last passage percolation, a change in the 
boundary densities a± can be compensated by a variation of the end point. For example, 
for the square [0,t>] 2 with boundary densities a± — (1 — v~ 1 ^ 3 w) ±1 , we can stretch the 1-axis 
by l/a+ and the 2-axis by l/a_. Then the boundary densities are 1 and the end point of the 
polymer is tilted by a distance v 2//3 away from (v, v), which is the quantity needed for the scaling 
form of the 2-point function. Thus for the PNG model the scaling form of K(6(x, t)0(O, 0)) 



follows from the analogue Proposition 4.1 for PNG and the asymptotic analysis in H, compare 



with Section 7.2 



The defect line along the diagonal is most easily visualized in the growth version of PNG, 
where unit up-steps are particles with velocity —1 and unit down-steps with velocity 1. If 
a_ = = a+, nucleation of pairs of steps are allowed only above the ground layer [—t, t]. Extra 
Poisson points along the diagonal with line density a correspond to a source at the origin which 
nucleates at rate a. The full PNG droplet has no reflection symmetry at the diagonal. With 
reflection symmetry one restricts oneself to the half droplet, where nucleations only above the 
ground layer [0, t] are allowed. In the depinned phase the extra mass is incorporated in the 
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droplet without visible modification of the macroscopic shape. In the pinned phase an extra 
pile of linear slope is created on top of the PNG droplet. 



6 Scaling theory 



In limit theorems one has to separate universal from model-dependent features. To give the 
standard example: let £j, j £ Z, be a stationary mean zero sequence of random variables which 
satisfies the central limit theorem as 



lim P 

n^oo 

" " 3=1 



Here Fq is the distribution function of a standard Gaussian random variable. Fq is universal 
(within well-studied limits), whereas the variance a 2 depends on the particular probability law 
for the £j's. In our example a 2 is easily computed as 

oo 

a 2 = E fe£o) • (6.2) 

j=-oo 

For the ASEP and similar growth models a corresponding scaling theory is available which 
determines the model-dependent longitudinal and transverse scales ||19|| . To apply the theory 



one has to know the stationary measures as labeled by their mean density, p. The two basic 
quantities are 

j(p) = average current at density p, 

A(p) = size of the density fluctuations in the stationary measure. 

If the density is regarded as the slope of a height function as in ( |2.1| ), then A(p) is the roughness 
amplitude for j \— > h(j) in the stationary measure, i.e. A(p) = lini^oo j _1 lE( [h(j) — h(0) — 
E{h(j)-h(0))} 2 ). 

For the TASEP j(p) = p p (r]o(l — r/i)) = p(l — p) and the height diffusion constant 

A(p) = lim -p p ((h(j) - (1 - 2p)j) 2 ) =4V(^( Wo ) - p 2 ) = 4p(l - p) . (6.3) 

3^<x> J 

For the discrete time TASEP, jump probability 1 — q, the stationary measure at density p 
is a Markov chain with transition probability Q, a 2 x 2-matrix which we write in the form 
Qoo = 1 — r, Q i = r, Q 10 = s, Qu = 1 — s. r,s are determined through (1 — r)(l — s) = qrs 
and (r + s)p = s. Then the average current 

j(p)=ps (6.4) 

and the roughness amplitude 

oo oo 

A(p) = 4 (E,( W )-p 2 )=4(p(l-p) + 2^(p(Q^) 11 -p 2 )). (6.5) 

j=-oo j=l 
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For the PNG model in the stationary measure the + , resp. — , particles are Poisson distributed 
with density p + , resp. p_, satisfying the stationarity constraint 2p + p_ = 1. The conserved 
density is p = p + — p_ and the current 

j(p)=P+ + P- = (2 + p 2 ) 1/2 . (6.6) 

From the Poisson property we conclude 

A(p)= P+ + p. = (2 + p 2 )^. (6.7) 

j'(p) = v(p) is the velocity of the density fluctuations and A = j"(p) is the KPZ coupling 
constant JIIJ. In terms of these quantities, the scale in the //.-direction is 

signA^IAIA 2 ) 1 /^ 1 / 3 (6.8) 

and the scale in the transverse j-direction is 

2 4/3 (A 2 A) l/3 t 2/3 ; (6g) 

For the TASEP, since A = —1, the /i-direction then comes in units of — 2 _1 / 3 (4p(l — p)) 2 / 3 t 1//3 , 
whereas the j-direction comes in units of 4(p(l — p)) 1 / 3 t 2 / 3 . 

In (|6.8|) and (|6.9| ) we have fixed two dimensionless scale factors. Their role is most easily 
explained in our entry example, where we could have adopted a definition of the error function 
Fq such that the Gaussian has mean 1/2. Then in (|6.1|) ax is to be replaced by y2ax indepen- 
dently of the value of a. The particular convention is determined through the comparison with 
one single test case. In our context we follow || in the definition of the distribution functions, 
who rely on the established conventions for the Painleve II Riemann-Hilbert problem. The test 
case ist the PNG model, which fixes the prefactors in ( |6.8|) to 2 -1 / 3 and in (|6.9|) to 2 4 / 3 . 

So far proofs are available only for a short list of models. In all cases the model-dependent 
parameters are such as to agree with the scaling theory. In particular the discrete time TASEP 
provides a two-parameter (p and q) test of the theory. 



7 Current fluctuations 

There is little doubt that the results of Baik and Rains also hold in the limit of an exponential 
distribution when q — > 1. Strictly speaking our results are conjectures, except for p_ = 1, 



p + = 0. We discuss the list of our examples. Some of them have been announced in 
where also numerical plots of F GUE , -^goe, and F are given. 



7.1 Initial step density 

Fluctuation results for N t are most easily summarized in the diagram of Figure |l|. We first 
remind the reader of the limit measure p^ for {rjtj, \j\ < iV}, iV arbitrary, t —>■ oo 0. The 
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GOE' 



P_ 




Figure 1: Fluctuation phase diagram for the TASEP. 



upper left corner is the maximal current phase with /i^ = pi/2- In the lower left p^ = p p _ and 
in the upper right p^ = Pp + - Along the line p_ + p + = 1 the long time limit is the superposition 
Poo = |pp_ + 2^p+- ^ n particular, we have E( N t ) ft —> ^00(^0(1 — Vi)) as t —* 00. 

To understand the fluctuations of N t we use the directed polymer picture and recall that 
the polymer starts at (0, 0) and ends at (n,n). If p_ < 1 — p+, p_ < |, then the polymer stays 
a finite fraction of steps at the right edge. Since for the edge steps the passage time is a sum 
of independent exponentials, its fluctuations are of order \/t and Gaussian (G). At some point 
the polymer must enter the bulk and we expect sub leading corrections of order t 1//3 , 

N t = p(l - p)t + t 1 ' 2 ^ - & 3 Z G jje (7.1) 

up to prefactors, with either p = p_ or p = p + . If p_ + p + = 1, p_ < |, then the polymer 
has a choice between the left and right edge and the fluctuations of N t are the maximum of 
two independent Gaussians (G 2 ). Recall that if £j has distribution function Fi, i — 1, 2, then 
their maximum has the distribution function FiF 2 provided £1,^2 are independent, hence our 
notation. 

If p_ > |, p + < ~, then it does not pay for the polymer to stay at neither edge. The 
fluctuations must be as for the Johansson case of trivial edges, p_ = 1, p + = 0, which have 
-^gue as distribution function. Thus the difficult cases are the two critical lines p_ = |, p + < |, 
resp. P- > |, p+ = |, and the critical point p_ = | = p + , which corresponds to the Bernoulli | 
initial measure. On the line p_ = | (p + = |) the directed polymer stays order n 1 / 3 at the left 
(right) edge, whereas at the critical point it makes a choice between the two edges. According 
to H the fluctuations of N t are of order t 1//3 and are GOE 2 distributed on the critical lines, F 
distributed at the critical point. 

From the point of view of the last passage percolation the end point (n, n) is somewhat 
special and one might consider more generally the end point (m,n) with the slope 9 = n/m 
fixed. To have a firm link to the TASEP we will use the parameterization Q2.1|) , in which case 
the quantity of interest is h t ([yt]), \y\ < 1, [•] denoting the integer part. h t (0) = 2N t is the 
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particular case just explained. The hydro dynamic theory [20, |T|] establishes the law of large 
numbers for h t ([yt]) with the result 

lim -M[yt\) = h(y) (7.2) 

t^co t 

almost surely for \y\ < 1. The limit h depends on p~,p+ and is given as follows: If p_ < p + , 
then 

- h{ )= f (l-2p_)y + 2p_(l-p_) for y<y e , 
\ (l-2p+)y + 2p + (l-p + ) for y > y c 

with y c = (p + (l - p + ) - p_(l - p-))/(p+ - p-). If p_ > p+, then 

(l-2p_)y + 2p_(l-p_) for y<l-2p_, 
%)=<! |(l/ 2 + l) for l-2p_ <y< l-2p + , (7.4) 

(1 -2p+)y + 2p + (l -p+) for l-2p+<y. 

We are interested in the fluctuations of 

h t ([yt])-fh(y) (7.5) 

for fixed y in the limit t — > oo. If ft, is linear around y, then the fluctuations are order \/t and 
Gaussian. If h has a cusp at y (and is linear on both sides), then the fluctuations are order \/t 
and G 2 . On the other hand, if h has nonzero curvature at y then the fluctuations are order t 1 ^ 3 
and GUE. The critical lines, GOE 2 , correspond to a y where at h(y) the curved piece joins the 
linear piece. Finally the critical point, F , is the merger of the two critical lines. 

We summarize our findings in the form of a (well-founded) conjecture. Pp_, P+ refers to the 
TASEP with p p _,p + as starting measure. h t (j) is defined in (|2~l|), h in ( |7.3|) , (|7.4|). 



Conjecture 7.1 (G) Let either p_ < p + and y > y c or p_ > p + and y > 1 — 2p + . Then 

lim ¥ p _ tP+ (th(y) - h t ([yt}) < (4p + (l - p + )(y - 1 + 2p + )t) 1,2 x) = F G (x) . (7.6) 

Let either p_ < p + and y < y c or p_ > p + and y < 1 — 2p_ . Then 

limF p _, p+ (th(y) - ht([yt}) < (4p_(l - p_)(-y + 1 - 2p_)t) 1/2 x) = F G (x) . (7.7) 

(G 2 ) Lei p_ < p + and y = y c; then 

lim P p _, p+ (t%) - Mfe/t]) < " P-)*) V V 

= F G {(4p + (l -p + ))- 1/2 x)F G ((4p_(l -p-))- l ' 2 x) . (7.8) 

(GUE) Let p_ > p + and 1 - 2p_ < y < 1 - 2p + . TTien 

lim P p _, p+ (f%) - /i t ([y*]) < 2~ 1 / 3 (1 - y 2 ) 2 / 3 * 1 ^) = F GUE (x) . (7.9) 

£—►00 \ / 
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(GOE 2 ) Let p_ > p + and either y = 1 — 2p_ or y = 1 — 2p + . T/ien 



lim P„ 



t^oo 



„ P+ (th(y) - h t ([yt]) < 2-V3(i _ ^2)2/3^/3^ 



(7.10) 



(F ) Let p_ = p = p + and y = 1 — 2p. Then 



lim P 



p- ,p+ 



th(y) - h t ([yt}) < 2-^(1 - y^'H^x] = F (x) . 



(7-11) 



i<G t/ie distribution function of a standard normal distributed random variable. F and 
the Tracy-Widom distribution functions -Fgue, -^goe are defined in J^/. 



The Gaussian case (G) with p_ = p + is proved in |2T| and the case (GUE) with p_ = 1, p + = 



m 



In ( |7.6|) to (|7.11| ) the random variable — (h t ([yt}) — th(y)) appears because of the inversion 
from the passage time G(m, n) to the height. In each case we have subtracted the asymptotic 
mean as computed from the hydrodynamic theory. For the variances of the Gaussians we have 
to determine how many steps, on the scale n, the directed polymer stays at either edge. For the 
prefactor 2 -1 / 3 (l —y 2 ) 2 ^ 3 in (\l.9\) to (|7. 1 1|) we either take the limit of the discrete time TASEP 



or use the scaling theory of Section |6|. Because of the inversion, also the lower and upper tails 
are interchanged. For example Fque(^) has the lower tail exp(— |x| 3 /12) and the upper tail 
exp(-4|x| 3 / 2 /3). 



7.2 Stationary two-point function 

The starting measure is Bernoulli p and P p is the corresponding path measure of the stationary 
TASEP, < p < 1. Density fluctuations propagate with velocity 1 — 2p. Therefore the 
anomalous fluctuations appear in ht([(l — 2p)t]), compare with ( |7.11|) . For the full two-point 
function we need the height distribution a distance of order t 2//3 away, i.e. at (1 — 2p)t + 4(p(l — 
p)) 1//3 t 2//3 w for arbitrary w, compare with ( |6.9| ). 

Conjecture 7.2 We have 

lim P p ( - h t ([(l - 2p)t + 4(p(l - p)) l 'H 2 l 3 w]) + 2p(l - p)t 

+ (1 - 2p) ((1 - 2p)t + 4(p(l - p^/H^w) < 2(p(l - p)) 2 / 3 ^) 
= F w (x). (7.12) 

The distribution function F w (x) is defined in j^J and denoted there by H(4w 2 + x; w, — w). The 
distribution function F appears also in ( \7.11\ ). 



Note that J dF w (x)x = 0, since we have subtracted the average of h t . The pref actors of x and 
w are infered from the scaling theory of Section |6|. For the PNG model, the analogue of ( |7.12j ) 
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Figure 2: A 3d-plot of the distribution functions F w for w > 0. 



follows from |J. If the convergence (|7.12| ) holds also for the second moments, then 

E,(fc([(l - 2p))t + 4(p(l - P )) 1/3 t^w]) - 2p(l - p)t 

-(1 - 2p)((l - 2p)t + 4(p(l - p^'Ww)] 2 ) 

^ 4(p(l - p)) 4 / 3 t 2 / 3 y dF w (x)x 2 

for large £. Setting 



and taking second derivatives in (|7.13|) , we conclude the scaling form of 

1.2 




1.2 1.4 



(7.13) 



(7.14) 



Figure 3: The second derivative of the second moment of F w vs. w. 

F w (x) is the solution of a first order partial differential equation with suitable boundary 
conditions. A numerical plot is provided in Figure 0. In essence it shows how the distribution 
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broadens for w — > oo, reflecting the cross over to the Gaussian fluctuations as one moves away 
from the characteristic {j = (1 — 2p)t}. The second moment, g(w), is symmetric and increases 
linearly as 4\w\ for large \w\. In Figure 3 we plot the scaling function g" as determined from 
the numerical evaluation of the second moment of F w at various values of w. The oscillations 
at small w result from numerical imprecisions for F w . In ]TI| an approximate nonlinear integral 
equation for g" has been derived, which was then solved numerically in [p^j . The approximate 
g(0) differs from the exact one by an order of 5%. From the available data a more accurate 
comparison does not seem to be feasible. 



7.3 Semi-infinite system 

We restrict to the half-lattice N and insert with rate a particles at site 1, respecting the 
constraint rji = 0,1. The initial measure is Bernoulli p p . The last passage percolation rep- 
resentation lives in the lower triangle {(i,j)\0 < j < i}. Along the diagonal the w(j,j) are 
exponential with rate a and at the lower edge the w(j, 0) are exponential with rate 1 — p. By 
maximizing the passage time on scale n we obtain the same phase diagram as for the critical 
step. However the fluctuations in N t , now the number of particles injected up to time t, are 



modified. In Figure 4 we summarize the findings |15| which could be written more formally as 
in Conjecture |7.1| and so far have been proven only for the Poisson case, i.e. the semi-infinite 
PNG model with a source at the origin. GSE is the distribution of the largest eigenvalue of a 
symplectic Gaussian random matrix. _p^ emi is a novel distribution and given by 

F scmi (x) = (1 + l(-v(x) + u{x)){x + 2u\x) + 2u{xf))F{x)E{xf (7.15) 

in the notation of |J. On the anti-diagonal a new feature appears. Recall that for the infinite 
system the distribution G 2 arose from the choice between the left and right edge. Here the 
directed polymer may stay for an arbitrary length at the lower edge and then switch to the 
diagonal. Optimizing over this one parameter family of paths yields for the passage time a 




Figure 4: Fluctuation phase diagram for the semi-infinite TASEP. 
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distribution which is given by the largest eigen value of a 2 x 2 GUE matrix, denoted by G 2x2 
in Figure [|. 

7.4 Defect lines, depinning 

For a change let us consider the Poisson last passage percolation. As already explained, in 
the square [0,v] 2 we have Poisson points with density 1. The directed polymer starts at (0,0) 
and ends at (v,v). We add extra Poisson points along the diagonal with line density a. Of 
primary interest is the last passage time, lim^oo v^Gaiv) = r(a). We know that r(0) = 2. 
The critical a c is defined as the largest a such that r(a) = 2. If a > a c , the polymer is pinned 
and typically returns after a length £\\(a) to the diagonal. Between two returns the polymer 
makes an excursion of size £±(a) away from the diagonal. The depinning transition is governed 
by r(a), £\\(a), £±(a) as a J, a c . If £±(a c ) < oo, the transition is first order, while £±_{oi) — > oo 
as a I a c signals as second order transition. 

Since there is considerable interest from the physics of disordered systems, a phenomeno- 
logical theory has been developed, which after some controversy seems to be widely accepted 
p3| , |24f . We explain the predictions and first recall the pure problem with the path measure 

(Z^P^exp [-P [ dsV(x s )] . (7.16) 
J o 

Here x s e l d is a Brownian path which has the a priori weight given by the Brownian bridge 
F BB with start and end point 0. V is a pinning potential, which by a suitable choice of units 
is simply V(x) = for \x\ > 1 and V(x) = —1 for \x\ < 1. Z t is the normalizing partition 
function. The analogue of r(a) is the ground state energy e(/3) = lim^oo — t^ 1 logZ t , e(0) = 0. 
For d — 3 a certain minimal strength is required to pin the Brownian motion. Thus (3 C > 
and e(/3) ~ — {(3 — /3 C ) 2 for (3 > f3 c . However for d = 2, Brownian motion is null-recurrent and 
j3 c = 0. Just a tiny bit of attraction suffices to pin the polymer. Correspondingly e(f3) has an 
essential singularity as e(f3) ~ — e -1 ^, £\\{/3) ~ e 1///3 , £±(/3) ~ £\\(j3) 1 ^ 2 for small j3. 

The claim is that a directed polymer in 1 + 1 dimensions in a random potential is pinned 
in essence as a directed polymer in 2 + 1 dimensions with zero bulk potential. Specifically, any 
attraction along the diagonal will pin the polymer, i.e. a c = 0, the passage time r(a) ~ 2+e _1//Q , 
and the longitidonal excursions £\\(q) ~ e 1//a with the usual link £j_(a) ~ £\\{a) 2 ^ for small a. 
No rigorous result seems to be available. 

Transcribed to the slow bond problem, the prediction is r c = 1 with an essential singularity 
for joo(r) at r = 1. 

To model the case where the Poisson points are reflection symmetric relative to the diagonal, 



we add in ( |7.16| ) a hard wall and replace V by Vh w , Vh w (xi, x±) = V(x) for x\ > 0, Vh w (x) = oo 
for Xi < 0. The entropic repulsion shifts to (3 C > even for d = 1,2. Therefore the prediction 
of the phenomenological theory is r(a) ~ 2 + c(a — a c ) 2 , £\\(a) ~ (a — ct c )~ 2 , and £±(a) ~ 
£||(a) 2 / 3 for a close to a c , a > a c . This is in agreement with the exact result a c = 1/a/2 and 
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r(a) = 2 for < a < a c , r(a) = \pl[a + (l/2a)) for a c > a, which slightly above a c yields 
r(a) ~2 + 2(a-a c ) 2 |3|. 

9 Conclusions 

It is rather surprising to have now, after more than 20 years of study, a technique available 
through which universal scaling functions can be computed, in some cases even very explicitly 
The main lesson to be drawn is that the universal distribution functions on the scale t 1 ^ 3 depend 
on the type of initial conditions. It would be of interest to understand whether our list is already 
complete. 

Despite progress there are several obvious elements missing. Firstly for the transverse 
fluctuations of the polymer we do not have such a fine information as on the passage time (= 
energy of the polymer) JTDJ. More importantly, only single distributions, like the passage time 
with given end points, can be handled. To have a more detailed understanding of the energy 
landscape joint distributions must be studied, like the joint distribution of the first passage 
times G(mi, rii), G(m 2 , n 2 ) refering to two distinct end points. 
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